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Abstract 

We make a precise scaling conjecture, based on renormalization group ideas, 

regarding the screening cloud around an impurity spin in the Kondo effect and 

test it numerically using the Density Matrix Renormalization Group method. 
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The renormalization group (RG) theory of the Kondo effect QTJ] suggests the existence of 
an exponentially large length scale, = Vp/T^ ~ ae l l pJ where vp is the Fermi velocity, 
Tk the Kondo temperature, a the unit cell dimension, J the Kondo coupling and p the 
density of states (per spin). One imagines a cloud of electrons of this size, typically microns, 
which screens the impurity spin. The meaning, or even the existence of this screening 
cloud has been a subject of some controversy [@-§J with no clear consensus emerging either 
experimentally or theoretically. In particular, the experiments of Boyce and Slichter were 
interpreted as indicating that the screening cloud had a size of order a lattice spacing, about 
a hundred times smaller than predicted by the renormalization group theory. In this letter 
we make a more precise statement of what the existence of this screening cloud really means, 
based on renormalization group scaling ideas, and test it against results from perturbation 
theory, local Fermi liquid theory and numerical simulations. A detailed account of our results 
will be presented elsewhere ||. Related theoretical work includes perturbative calculations 
HQ, RG approaches and numerical work on the Anderson model ||. For a general 
review of the Kondo effect see Ref. [[| and references therein. 

We consider the standard Kondo model: 

H = £ ^ k >ka + JS imp ■ £ Wc^T/W, (1) 
k k,k' A 

with Si mp = 1/2. We are interested in the Knight shift, proportional to the local suscepti- 
bility: 

X (r,T) = (1/T) < ^WyVW^ot > . (2) 

Here Stot is the total spin operator, including both impurity and conduction electron spins; 
we set the two g-factors equal. In the scaling limit, rkp ^> 1, T <C Ep, we conjecture the 
scaling form: 

X (r,T, J) - t = -L, {/ (H, p) cos Ur + S (H, £)1 +g (H, £)} , (3) 

where /, g and <5 are universal functions of two scaling variables, and the bulk Pauli suscep- 
tibility has been subtracted. The prefactor of 1/r 2 arises from reducing the 3 dimensional 



problem to an effective one- dimensional problem, |I| which is assumed to obey scaling. 



Instead of the second scaling variable, T/Tk, we may equivalently use the renormalized 
Kondo coupling at scale T, A e ff(T), where A is the dimensionless Kondo coupling, A = pj. 
For weak bare coupling and high temperatures, T/Tk ~ exp[l/A e fr(T)]. Eq. (|3|) follows 
naturally from the expected asymptotic limits of \ and standard scaling hypotheses. This 
scaling form is probably the best definition of what it means to have a screening cloud since 
it implies that for T < Tk, x{ r ) varies over at least a distance of £k- 

Eq. (|) is consistent with perturbation theory at T ^> Tk, where X e «(T) is small. The 
first order result is: 

v(r, T) — ~ = „ rp cos2k F r, (4) 

AV ' 2 16r 2 w F sinh^ ' v ; 

consistent with Eq. (|]). Including the next order term, |6]||] and assuming also, rT/vp <^ 1 
we find: 

X(r) ~ % = - C " S2 ^ r { A + A 2 [ln(fc F r) + constant]}. (5) 

In low order perturbation theory, the effective coupling at temperature T is given by |]J 

Aeff(T) = A + A 2 ln(^ F /T) + 0(A 3 ). (6) 

Thus, to 0(A 2 ), the quantity in brackets in Eq. (|D may be written, as 

{A cfr (T) + X cS (T) 2 [\n(rT/v F ) + constant]}, (7) 

consistent with Eq. (|3|). 

Conversely, at T Tk and r ^> £k, we expect local Fermi liquid theory |J to apply. 
X then reduces to the magnetic susceptibility of a non- magnetic impurity with a 7r/2 phase 
shift. This can be obtained by differentiating the Friedel oscillation formula for the local 
charge density: 

n(r) = n - \ 3 cos[2k F r + 7r/2], (8) 
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with respect to the chemical potential yielding: 

/ ^ 1 dn p 1 , , 

again consistent with Eq. (|3p. Note that x{ r i T) is longer-range at low T after the screening 
cloud has formed, \ oc 1/r 2 , than at higher T before it has formed, \ oc 1/r 3 , (Eq. (|5])). 
Integrating Eq. @ over d 3 r, to obtain the long-range contribution to the total susceptibility, 
only g contributes, giving (after a change of variables): 

Air f°° 

Xlong range = 777 / dxg(x,T/T K ). (10) 
1 JO 

This is consistent with the low temperature Fermi liquid prediction |IJ x ~ 1/Tk provided 
that / °° dxg(x, y) oc y (for y <C 1). 

The functions f(x,y), g(x,y), and 8(x,y), in the regime y < 1, x/y = t/^x < 1 are 
of special interest. They describe the interior of the screening cloud at low T. One might 
naively suppose that a small r <C £^ would also cut off the renormalization of the effective 
coupling so that deep inside the screening cloud we recover weak-coupling behavior (for weak 
bare coupling) even at low T. However, explicit RG calculations by Gan ||[] show this not to 
be the case. He showed that a small r does not cut off the infrared divergence of the Kondo 
coupling in perturbation theory. Thus the scaling functions are expected to be non-trivial 
in this region. 

Now let us consider the experiments of Boyce and Slichter on Fe doped Cu. They 
measured the Knight shift from different shells of Cu atoms within a few lattice spacings 
from the impurities from T=300K down to well below the Kondo temperature of 29K. They 
found the factorized form: x(r, T) = f(r)/{T + T^),with some rapidly varying function 
fir) (which, in fact, changes sign over the small range of r considered). Note that all 
measurements are taken in the regime rT/vp -C 1, rT^/vp <C 1. In fact the values of r are 
so small that it is unclear whether the scaling form of Eq. ([$]) holds at all. Assuming it does, 
we may essentially consider x — > in the scaling functions. Assuming that, in this limit, / 
dominates over g, and that S is essentially constant, we can obtain the above result if f(x, y) 
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exhibits an approximately factorized form for x/y <C 1 and all y: f(x,y) ~ f{x/y)/(y + 1). 
Note that such factorization could not also occur at large r, r ^> if our assumed scaling 
and asymptotic behaviors are correct. In this region, the 1/r 3 behavior of the 2k f part at 
T ^> Tk crosses over to 1/r 2 at T Tk- Thus, a conclusive test of the existence of the 
screening cloud would probably require experiments which probe length scales of O 

To test our scaling hypothesis using the density matrix renormalization group (DMRG) 



method [JIT] we consider a tight-binding model with open boundary conditions: 



with the total number of sites, L, even. We shall always take the hopping strength, t, to be 
unity, and we only consider ground-state properties at half-filling. We keep m = 128 states 



and for optimal precision we apply the finite lattice DMRG JXTJ] to each step in the infinite 
lattice DMGR method. In the ground-state the total magnetization, S* ot = 1/2, is a constant 
and we therefore consider the local magnetization, < >, instead of A% = x( r , T, J) —p/2. 
Here is the z-component of the electron spin at site j. With the DMRG method we can 
obtain results for finite systems, L < 40 — 50, at T = 0. Hence we need a finite-size form 
of Eq. (JJ). This can be obtained by noting that since kp = n/2 at half-filling, 5 can be 
eliminated from Eq. ([|). Substituting rT '/vp — > r/L and T/Tk — > £k/L we get: 

«*>-HKv& l -» + 'H)Y (12) 

In writing this equation we have included a factor of T/r 2 in going from 3D to ID and 
converting to the local magnetization. 

We briefly mention a few useful results regarding the weak and strong coupling behavior 
of Eq. (JIl]) . In the weak coupling limit, J <C t, using perturbation theory, we find for large 
L: 

< s* >^ ^-(-iy. (13) 

In the strong coupling limit, J ^> t, it can be shown that the on-site magnetization behaves 
in the following way: 



<^>=I s in 2 ^, j>2. (14) 

These results are consistent with the scaling form Eq. ([X2|) if we note that in the weak 
coupling limit we can substitute J/t with the effective coupling l/ln(£^/L). 

We start by considering < / 2 >. Restricting our analysis to even values of L/2 Eq. ( [3~2|) 
takes the simpler form L < S Z L , 2 >= h(L/^x)- If this scaling form is obeyed it should be 
possible to collapse data for different L and J onto a single scaling curve. In Fig. [I] we show 
our results for a large range of Kondo couplings J with £jc(J = 2.5) fixed to 1. Choosing 
appropriate values of all the data can clearly be brought onto the predicted scaling form. 
A very dramatic crossover from weak to strong coupling behavior is evident. Note that for 
strong coupling L < S^^ > approaches 1 consistent with the result Eq. (|T4]) . The weak 
coupling result Eq. (|T3"D predicts a roughly constant value of < Sf j , 2 > for fixed J consistent 
with the behavior in Fig. [j]. The fact that < 5| >— > constant/L for j 3> £k is analogous to 
the prediction that local susceptibility, x(r, T), should be longer range when the screening 
cloud has formed than when it has not, as discussed above. 

We proceed to analyze the on-site magnetization for fixed L/£k as a function of j. 
Our results are shown in Fig. (Q) where we display results for J = 1.5 and J = 1.8 with 
£ K {J = 1.5) ~ 4.85 and £ K (J = 1.8) ~ 2.7, and we use L = 36,20 for J = 1.5,1.8 
respectively. Thus, L/^k ~ 7.4 is kept fixed. An almost complete data collapse is clearly 
visible. This is a highly non-trivial test of the scaling relation Eq. (|T2|) and hence also of 
Eq. (JJ). Analogous scaling plots can be performed for the other couplings studied. 

Having extracted £x from the scaling of L < S z L j 2 >, we can study as a function of 
J . As mentioned in the introduction, usual renormalization group arguments predicts to 
lowest order an exponential form for Including one more term in the renormalization 



group equation for the effective coupling we have fT2|j . 

-dA/dlnA = A 2 -A 3 /2. (15) 

Here A is the momentum space cut-off. We can integrate this equation to get an improved 
expression for This gives 
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6r = foe 



l/Ao-l/c 



\ (16) 



Here c is positive and 0(1). The inverse of the bare cut-off, £o, should satisfy 1 <C £o £k 
and Ao, the bare coupling, is given by Ao = Jp+(Jp) 2 ln[tan(2£ /a)] with the density of states 
for the tight-binding model, p = l/2irt. This form is valid for small Ao- The corrections 
arising from the constant c is presumably only one of several equally important terms; we 
have included it here to improve the agreement with the numerical results. In the limit 
A -> Eq. ([ID reduces to: ( K oc e 1/x °/y/\^, in agreement with previous results ||13|| . The 



form, Eq. ( |T6"|) is obtained using weak coupling RG and should therefore only be valid for 
the range of J where 1 C C L. In this region we can now try to fit our previously 
obtained results for to the form Eq. fllTf). This is shown in Fig. |3| where the DMRG 
results are shown as circles and the solid line indicates a least square fit to Eq. ([If]) with 
the parameters £o = 1.76, c = 0.21. For these intermediate couplings the RG form works 
surprisingly well, and the data seems clearly consistent with the expected form. However, 
one should bear in mind that for weaker couplings £k has sizable error bars that we cannot 
estimate; secondly, the extracted £k is only determined up to a multiplicative constant. 

In conclusion we have presented several non-trivial test showing that the scaling form 
Eq. (|H) clearly is obeyed in its finite-size form. Finite-size effects may be important in recent 
experiments using samples with lengths of order £^ or smaller ||. The associated Kondo 
length scale, which roughly determines the size of the screening cloud diverges at weaker 
couplings with a behavior consistent with higher order renormalization group calculations, 
Eq. ([TED . The Knight shift associated with the local susceptibility is longer range at low 
temperatures where the screening cloud has formed than at higher temperatures where it 
has not. Inside this screening cloud, as shown by Gan Q], one does not expect weak coupling 
behavior to be recovered. 
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FIGURES 
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FIG. 1. L < S Z L , 2 > as a function of L/^k for a range of different coupling constants. The 
initial point corresponds in all cases to L = 4, subsequent points have L increased by 4. 
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FIG. 2. L times the expectation value of the z-component of the electron spin, < S? >, 
as a function of ]/£,k{J)- Two systems are shown: J = 1.8, £x = 2.7, L = 20 and 
J = 1.5, = 4.85, L = 36. Thus in both cases we have L/^k ~ 7.4. 
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FIG. 3. The Kondo length, as a function of Kondo coupling, J. The circles denote the 
numerical DMRG results. The solid line indicates a least square fit of ^k(J) for J = 1 — 4 to the 
form Eq. ([IB]). The fitted parameters are £o = 1.76, c = 0.21. 
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